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LONG-TERM  GOALS 

The  long-term  goals  are  to  contribute  improvements  in  eurrent  physieal  understanding  and  modeling  of 
interfaeial  proeesses  fundamental  to  air-sea  interaetion  fluxes,  partieularly  those  involving  wave 
breaking  and  spray  droplet  produetion.  These  advanees  will  improve  the  reliability  of  operational  sea 
state  and  oeean  weather  foreeasting  models,  partieularly  for  severe  to  extreme  sea  states. 

OBJECTIVES 

This  projeet  seeks  to  improve  the  reliability  of  air-sea  interfaeial  flux  parameterizations  in  coupled  sea 
state/marine  weather  foreeasting  models,  with  a  particular  focus  on  refining  and  incorporating  the  role 
of  wave  breaking  and  sea  spray  in  severe  eonditions.  Our  approaeh  is  through  developing  more 
realistie  parameterizations  for  breaking  oeeurrenee  and  strength,  and  sea  spray/spume  souree  funetions 
and  validating  them  in  our  test-bed  models.  Our  goal  is  then  to  implement  these  improvements  in  a 
coupled  COAMPS/WaveWateh  III  model  for  operational  use. 

APPROACH 

Our  approaeh  is  to  build  substantially  on  our  aceumulated  expertise  in  sea  surface  proeesses  and  air-sea 
interaetion  (Banner)  and  numerieal  weather  modelling  (Leslie)  to  identify  and  fill  fundamental 
knowledge  gaps  in  order  to  improve  the  modelling  aeerrraey  for  severe  marine  meteorologieal  events 
sueh  as  hurricanes.  Such  improvements  depend  eritically  on  gaining  a  more  eomplete  understanding  of 
severe  sea  state  phenomena  linked  to  wave  breaking  -  inereased  surfaee  drag  and  sea  spray  production. 
Quantifying  the  distribution  of  wave  breaking  events  has  progressed  within  this  projeet  as  a  result  of  an 
intensive  eollaboration  with  D.  Farmer  and  J.  Gemmrieh  (lOS,  Canada)  in  FYOl.  From  their  storm 
wave  datasets,  a  robust  threshold  behavior  was  identified  in  terms  of  the  local  wave  spectral  saturation 
level  B(/)  for  wave  breaking  at  different  scales  (Banner  et  al,  2002,  hereafter  BGF02).  The  speetral 
saturation  B(/)  =  fY{f)l2g  for  the  wave  energy  spectrum  Y(f).  This  result  has  proven  very  useful 
in  our  ongoing  effort  to  parameterize  wave  breaking  speetrally  in  wind-wave  models,  and  underpins 
our  new  spectral  source  term  formulations  for  the  wave  energy  dissipation  rate.  It  is  also  being  used  in 
the  spray/spume  production  term  we  are  formulating  in  eollaboration  with  C.  Fairall  (NOAA). 
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Our  model  development  has  foeused  on  refining  full-bandwidth  eomputations  of  the  direetional  wave 
speetrum  and  its  tail  region  using  ‘exaet’  versions  of  the  nonlinear  wave-wave  interaetion  souree  term 
in  the  radiative  transfer  equation  for  the  wave  field.  We  have  coneentrated  on  feteh-limited  growth 
oases,  as  this  is  where  data  exists.  Suoh  an  approaoh  is  necessary  to  ensure  that  the  modelled  spectral 
saturation  levels  and  dissipation  rates  are  consistent  with  observed  levels.  This  is  essential  for 
predicting  spectral  breaking  wave  properties,  and  for  reliable  calculations  of  the  enhanced  fluxes 
between  the  atmosphere  and  ocean  associated  with  wave  breaking  events.  A  successful  outcome  in 
this  context  is  needed  to  underpin  future  operational  versions  of  our  spectral  wave  breaking  models. 
Our  effort  has  focussed  strongly  on  formulating,  implementing  and  refining  modeling  strategies  for  (a) 
extracting  the  relevant  wave  breaking  parameters,  (b)  calculation  of  wind  stress/roughness  length 
enhancements  and  updating  the  surface  layer  winds  accordingly.  For  (a),  we  have  significantly  refined 
our  capabilities  reported  previously  for  calculating  the  spectral  density  of  mean  breaking  wave  crest 
length/unit  area.  This  is  one  of  the  primary  goals  in  our  modeling  effort,  as  this  quantity  is  central  to 
the  prediction  of  breaking  wave  enhancements  to  the  wind  stress,  and  development  of  a  spray/spume 
source  function  based  on  sea  state  rather  than  the  wind  field.  Our  progress  builds  on  refinement  of  the 
form  of  Sds  proposed  by  Alves  and  Banner  (2002)  [hereafter  AB02]  and  on  advances  in  predicting 
wave  breaking  probabilities  at  different  wave  scales  in  BGF02.  The  AB02  form  of  Sds  was  upgraded 
to  incorporate  the  observed  BGF02  breaking  saturation  threshold.  We  have  also  investigated  various 
refinements  to  its  spectral  distribution  in  order  to  provide  a  better  match  to  the  wind  input  source 
function  Sin  at  higher  wavenumbers,  and  to  recently  available  published  spectral  wave  breaking 
observations  of  Melville  and  Matusov  (2002)  [hereafter  MM02].  To  progress,  we  have  had  to  invest 
very  considerable  effort  evaluating  various  versions  of  the  ‘exact’  Sni  code  and  propagation  schemes  to 
ensure  accuracy  and  minimise  computational  instabilities  that  can  develop  at  higher  wavenumbers. 

WORK  COMPLETED 

1 .  Wind- wave  evolution  modelling 

In  brief,  during  FY03  we  continued  our  main  effort  on  model  development.  We  refined  our  capability 
for  computing  the  directional  wave  spectrum,  with  particular  emphasis  on  the  spectral  tail  region  out  to 
frequencies  of  up  to  3Hz.  This  involves  large  numbers  of  very  lengthy  computations,  associated  with 
the  well-known  complexity  of  the  exact  form  of  Sni.  We  concentrated  our  attention  on  wind  speeds  in 
two  ranges  (a)  7-14  m/s  and  (b)  30  m/s.  For  (a),  direct  comparisons  were  made  with  the  ID  transect 
wave  spectra  and  the  breaking  crest  length  spectra  of  MM02,  and  recent  results  are  presented  below. 
For  (b),  the  aim  was  to  demonstrate  stable  model  tail  behavior  for  gale  force  wind  speeds,  which  was 
achieved.  A  detailed  investigation  was  made  of  the  new  analysis  technique  (see  2.  below)  that  we 
developed  for  determining  wave  crest  passage  rates  from  ocean  wave  data  needed  for  our  modeling 
effort  on  breaking  waves. 

In  greater  detail,  we  formulated  a  number  of  variants  of  the  unknown  spectral  wave  energy  dissipation 
rate  source  term,  Sds,  each  with  a  nonlinear  dependence  on  the  saturation  ratio.  SdO  and  Sd2  are  two 
variants  with  different  wavenumber  weightings.  We  investigated  these  in  the  spectral  wave  evolution 
equation  for  fetch-limited  conditions,  using  a  full  spectral  bandwidth,  i.e.  free  of  any  imposed  spectral 
tail.  The  nonlinear  transfer  source  term  Sni  used  was  the  standard  Resio-Tracy  exact  form,  as  this  gave 
the  best  stability.  The  wind  input  source  term  Sm  was  based  on  a  modified  parameterisation  of  Yan 
(1987).  It  embodies  the  traditionally-  observed  dependences  on  wind  parameters  linear  near  the 
spectral  peak  and  quadratic  in  the  spectral  tail.  It  also  provides  input  to  the  wind  for  waves  outrunning 
the  wind.  Modeled  directional  spreading  properties  compare  very  favorably  with  the  data  of  Hwang  et 
al.  (2000). 


Throughout  these  ealeulations,  our  results  for  integral  wave  speetral  properties  sueh  as  mean  energy 
and  peak  frequeney  aehieve  a  very  elose  eorrespondenee  with  the  observed  evolution  behavior  in 
Kahma  and  Calkoen  (1992).  We  also  obtain  very  elose  model  agreement,  both  level  and  speetral 
slope,  with  the  high  wavenumber  tail  of  the  1-D  transeet  speetrum  observed  by  MM02  (not  shown 
here).  Following  on  from  these  validations,  we  made  extensive  ealeulations  of  the  speetral  density  of 
breaking  erest  length  per  unit  area,  and  how  these  breaking  wave  speetra  ehange  with  wind  speed  and 
wave  age,  whieh  are  partieularly  demanding  tests.  We  eompare  our  model  results  with  the  available 
observational  data  of  MM02,  with  very  mature  wind  sea  eonditions  and  winds  of  9.8  m/s  and  13.6  m/s, 
From  the  speetral  evolution  ealeulations,  of  partieular  importanee  are  the  speetral  distributions  of 
dissipation  rate,  Sds,  speetral  saturation,  B,  and  the  breaking  erest  length/unit  area  of  waves  travelling 
with  speeds  in  (e,  e+de),  A(e).  Our  ealeulations  highlight,  for  the  first  time,  the  sensitivity  of  the  A(e) 
predietions  to  the  imposed  Sds-  As  there  are  no  available  direet  speetral  dissipation  rate  observations, 
the  results  of  MM02  for  A(o)  provide  the  first  independent  measurement  that  ean  be  used  to  deeide 
between  Sds  formulations.  This  has  underpinned  our  systematie  investigation  of  the  impaet  of 
eompeting  Sds  formulations.  For  sueh  a  eomplex  nonlinear  system,  this  represents  a  very  ehallenging 
task.  Progress  on  this  fundamental  problem  is  reported  under  ‘Results’. 

2.  Riding  wave  analysis 

For  1.,  we  developed  a  novel  time  series  analysis  teehnique  (‘riding  wave  removal’  [RWR]  method) 
for  eounting  riding  wave  erest  passage  rates  needed  to  model  breaking  waves.  This  teehnique  eounts, 
and  also  eharaeterises  the  geometrieal  properties  (e.g.  steepness)  of,  the  short  wind  waves  riding  on 
longer  waves.  We  investigated  these  properties  for  a  number  of  field  data  sets  (North  Sea,  Blaek  Sea, 
N.  Paeifie)  for  young  and  old  wind  seas.  This  teehnique  ean  provide  useful  ‘physieal  roughness’  data 
eomplementary  to  the  usual  speetral  data,  and  a  paper  is  in  preparation. 

3 .  New  wind  input  souree  funetion 

A  new  wind  input  souree  funetion  has  been  proposed  (Donelan  et  al,  2003).  This  is  based  on  Lake 
George  data  for  very  young  sea  eonditions,  whieh  shows  an  abrupt  transition  of  the  near-surfaee  air 
flow  from  attaehed  to  separated,  when  the  waves  beeome  sufficiently  steep.  This  is  accompanied  by  a 
marked  reduction  in  the  input  momentum  and  energy  fluxes  from  the  wind.  However,  the  strongly 
discontinuous  nature  of  this  form  of  wind  input  term  presents  new  modeling  challenges  that  we  felt 
should  await  further  progress  with  1 .  before  implementation. 

4.  Improved  spume/spray  droplet  source  term 

In  collaboration  with  two  other  CBLAST  Pi’s:  C.  Fairall  and  W.  Asher,  an  exploratory  laboratory 
experiment  was  staged  at  the  UNSW  Water  Research  Laboratory  in  Sydney,  Australia.  The  primary 
aim  was  to  seek  a  dependence  of  spume  droplet  flux  on  breaking  wave  properties  rather  on  the  wind 
speed.  This  is  the  basis  of  an  initial  wave  breaking/turbulent  flow  model  of  spume  droplet  production 
in  terms  of  the  surface  dissipation  rate  of  breaking  waves  developed  with  C.  Fairall.  Droplet 
measurements  were  made  with  Fairall’s  CIP  probe  and  with  Asher’s  TSI  Droplet  analysis  system  as  a 
precursor  to  their  aircraft  deployment  in  the  CBLAST  Hurricane  missions.  Typically,  our  runs  had 
friction  velocities  in  the  range  1.2-1. 8  m/s  in  conjunction  with  paddle-initiated  wind  waves  with 
frequencies  in  the  range  1 .36-1 .6  Hz.  Water  salinity  was  varied  from  0  to  25  parts  per  thousand.  The 
breaking  wave  properties  were  measured  with  a  small  directional  array  of  surface-piercing  wire 


impedance  gages  each  sampled  at  over  2000  Hz,  in  order  to  capture  the  surface  fine  structure.  Overall, 
we  successfully  gathered  a  comprehensive  data  set  on  droplet  properties  and  the  associated  wave 
surface  microstructure.  Analysis  of  the  breaking  wave  data  to  extract  a  consistent  measure  of  ‘surface 
disturbance  energy  levels’  is  complex  and  currently  in  progress. 

5.  Coupling  wind  profile  and  sea  surface  aerodynamic  roughness  changes 

A  Fortran  code  has  been  developed  to  allow  the  iterative  computation  of  the  wind  stress  and  wave 
spectrum  in  response  to  changes  due  to  breaking-induced  wave  drag.  This  approach  is  based  on  the 
calculation  of  wave  breaking  properties  from  the  wave  spectrum,  the  consequent  enhanced  wave  drag 
and  the  adjustment  of  the  aerodynamic  roughness  length  zq  in  the  assumed  logarithmic  mean  velocity 
profile  for  the  surface  layer  wind  field.  In  turn,  the  wave  spectrum  is  modified  in  response  to  the 
updated  wind  profile.  This  code  awaits  introduction  into  the  coupled  test-bed  model  once  the  breaking 
wave  extraction  algorithm  in  the  model  has  been  optimised  against  available  data. 

RESULTS 


Wind  wave  modelling  and  spectral  breaking  wave  distributions 

Fig.  1  shows  the  excellent  correspondence  between  modeled  and  observed  evolution  for  the  key 
integral  wave  properties,  using  SdO  as  the  form  of  Sds-  Similar  results  were  obtained  with  Sd2.  Fig. 2 
shows  the  substantial  wind  speed  and  wave  age  variation  resulting  from  our  Sd2  dissipation  rate  source 
term  spectrum.  The  wind  speeds  match  those  in  the  old  wind  sea  observations  of  MM02. 


Fig.  1  Comparison  of  modeled  and  observed 
Kahma-Calkoen  (1992)  integral  wave  properties 
for  fetch-limited  evolution,  using  our  saturation 
threshold  dissipation  rate  source  term.  Non- 
dimensional  wave  energy  e*  (upper  panel)  and 
peak  frequency  v*  (lower  panel)  against  non- 
dimensional  fetch  based  on  u*  scaling. 


Fig.  2  These  panels  show  the  strong  wind 
speed  variation  of  modeled  dissipation  rate  Sds 
spectra  computed  using  Sd2  for  Uio  =  9.8m/s 
and  13.6  m/s  for  young  wind  seas  (upper  panel) 
and  old  wind  seas  (lower  panel).  The 
horizontal  axis  is  frequency  normalized  by  the 
frequency  at  the  peak  of  the  wave  spectrum. 


It  is  only  when  the  key  speetral  tail  properties  are  eompared  can  large  differences  be  seen  between  the 
Sds  source  terms.  Fig. 3  below  shows  the  spectral  saturation  B  for  [/io~10  m/s  using  the  two  forms  SdO 
and  Sd2,  as  a  function  of /%,  compared  with  the  saturation  data  from  the  North  Sea.  Sd2  captures  the 
levels  and  overall  shape  of  the  observations,  while  SdO  shows  a  non-physical  decrease  above  the 
spectral  peak.  Fig. 4  shows  the  sensitivity  to  the  choice  of  Sds  of  our  model  computations  of  A(c),  the 
spectral  density  of  breaking  crest  length/unit  area  of  sea  surface,  expressed  as  a  function  of  the  phase 
speed  c.  The  sea  state  conditions  matched  those  in  the  observations  of  MM02.  The  computed  results 
using  Sd2  follow  the  spectral  steepness  of  the  MM02  observations  more  closely.  We  note  that  MM02 
only  include  growing  breaking  crests  in  their  analysis,  and  hence  neglect  the  dissipation  associated 
with  relaxing  breaking  crests.  This  is  consistent  with  our  systematically  higher  spectral  levels 
compared  with  MM02.  A  key  issue  is  to  optimize  the  wave  age  and  wind  speed  performance  of  Sds  in 
predicting  A(c).  Sea  state  data  gathered  during  the  recent  CBLAST  hurricane  missions  will  be  crucial 
to  this  goal. 


Fig.  3  Comparison  of  computed  azimuthally- 
integrated  spectral  saturation  B  as  a  function 
of  f/fp  for  U10  =  9.8  m/s  (upper panel)  and 
13.6  m/s  (lower  panel)  using  two  versions  (SdO 
and  Sd2)  of  the  spectral  dissipation  source 
term  Sds.  A  typical  saturation  spectrum  for  9.8 
m/s  observed  in  the  North  Sea  is  shown  for 
comparison. 


Fig.  4  Computations  of  A(c),  the  spectral  density 
of  breaking  crest  length/unit  area  as  a  function 
of  phase  speed  c,  using  two  versions  (SdO  and 
Sd2)  of  the  spectral  dissipation  source  term  Sds- 
Wind  speeds  Uio  are  9.8  m/s  (upper panel)  and 
13.6  m/s  (lower  panel)  and  seas  are  very  old 
(Uio/Cp~0.8).  The  MM02  observations  for  each 
wind  speed  are  shown  for  comparison. 


IMPACT/APPLICATIONS 

Enhanced  scientific  understanding  of  severe  sea  state  air-sea  interfacial  processes,  particularly  wave 
breaking  and  spray/droplet  production  rates,  will  provide  more  reliable  parameterizations  of  these 
processes  and  closely  related  air-sea  fluxes.  These  improved  parameterizations  will  increase  the 
reliability  of  operational  sea  state  and  marine  meteorological  forecasts,  especially  during  severe  marine 
weather  conditions.  A  particular  benefit  for  smaller  vessels  will  be  routine  forecasts  of  occurrence  rate 
of  dangerous  breaking  waves. 


RELATED  PROJECTS 


The  ONR  project  Source  Term  Balance  for  Finite  Depth  Wind  Waves  (Young,  Banner  and  Donelan) 
includes  a  strong  focus  on  the  influence  of  steep  waves  and  wave  breaking  on  the  wind  input  source 
function  in  strongly  forced,  constant  depth,  shallow  water  environments.  This  data  set  has  been 
analyzed  and  initial  papers  on  the  methodology  and  results  have  been  submitted,  with  other  papers  in 
preparation  [e.g.  Donelan  et  al,  2003a,b;  Young  et  al,  2003].  We  have  begun  exploring  these  outcomes 
in  this  project. 
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